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Abstract 

We present a formulation for the construction of first order equations which 
describe particles with spin, in the context of a manifestly covariant relativistic the- 
ory governed by an invariant evolution parameter; one obtains a consistent quan- 
tized formalism dealing with off-shell particles with spin. Our basic requirement 
is that the second order equation in the theory is of the Schrodinger-Stueckelberg 
type, which exhibits features of both the Klein-Gordon and Schrodinger equations. 
This requirement restricts the structure of the first order equation, in particular. 
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to a chiral form. One thus obtains, in a natural way, a theory of chiral form for 
massive particles, which may contain both left and right chiralities, or just one of 
them. We observe that by iterating the first order system, we are able to obtain 
second order forms containing the transverse and longitudinal momentum relative 
to a time-like vector i^i'^ = — 1 used to maintain covariance of the theory. This 
time-like vector coincides with the one used by Horwitz, Piron, and Reuse to ob- 
tain an invariant positive definite space-time scalar product, which permits the 
construction of an induced representation for states of a particle with spin. We 
discuss the currents and continuity equations, and show that these equations of 
motion and their currents are closely related to the spin and convection parts of 
the Gordon decomposition of the Dirac current. The transverse and longitudinal 
aspects of the particle are complementary, and can be treated in a unified manner 
using a tensor product Hilbert space. Introducing the electromagnetic field we 
find an equation which gives rise to the correct gyromagnetic ratio, and is fully 
Hermitian under the proposed scalar product. Finally, we show that the original 
structure of Dirac's equation and its solutions is obtained in the highly constrained 
limit in which is proportional to on mass shell. The chiral nature of the the- 
ory is apparent. We define the discrete symmetries of the theory, and find that 
they are represented by states which are pure left or right handed. 
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I. Introduction 



The Dirac theory of the spin ^ particles, hke the Klein-Gordon description of particles 
without spin, makes use of wave functions which are not localized, as pointed out by 
Newton and Wigner The equations of motion, although useful for describing the 
properties of quantum fields, are not adequate for the construction of effective one- 
particle quantum theories. It has been found that the use of an invariant parameter to 
describe the evolution of states which are off-mass-shell provides a general framework in 
which this problem is solved [0, 

First order equations for particles with spin ^ , for which the evolution of the entire 
system is governed by an invariant parameter, have been searched for and studied by 
many over the years. Their importance lies in describing the behavior of the spin degrees 
of freedom of fermionic particles in such theories. We can trace the formulation of 
covariant quantum theories with an invariant evolution parameter from Fock [Q, through 
Stueckelberg [Q, Nambu 0, Schwinger and Feynman up to more recent work 
done by Cooke |^ , Horwitz and Piron ||^ , and Fanchi []TU| . 

Second order equations for spin | particles have been found and studied by Fock 
Feynman [Q, Horwitz et al |]TT|, and Reuse |T^. A number of first order equations for 
spin I particles have been proposed by Nambu , Feynman [|| , Kubo , and Davidon 
[p!4| . The Dirac equation [|l^ does not contain an invariant evolution parameter, and 
applies to a three-dimensional measure space. For a summary of the subject and an 



extensive list of references see Fanchi Some of the first order equations were intro- 
duced in an ad hoc manner, and some of the formulations were incompatible with the 
postulated second order equation for the evolution of a free particle in the corresponding 
theories. In some cases they led to free evolution equations which are second order in the 
invariant parameter, admitting solutions which propagate forward and backward, thus 
invalidating the interpretation of the invariant parameter as an unidirectional evolution 
parameter |T^, a feature that is quite important for the interpretation of the theory. 

In this paper we propose a first order equation of motion for a spin | particle, in 
the framework of the formalism developed by Horwitz and Piron , which is consistent 
with the form of the second order evolution equation for a free particle, the Schrodinger- 
Stueckelberg equation 

id,^ = (1) 
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It seems that this kind of evolution kernel, proportional to P^P^, is best suited to 
describe a covariant theory (see Fanchi |T9|). It leads to the correct relation for the 
velocities of a classical particle and permits separation of variables in the many body 
case [|]. 

We find that to achieve the goal of obtaining a first order equation for a spin | particle, 
it is necessary to introduce nilpotent operators which result in a chiral theory even though 
we are not restricted to massless fermions. This result may be relevant to the present sit- 
uation in weak interactions and neutrino physics, where one finds chiral fields although 
it is not clear the neutrino mass is zero, and when theories beyond the standard model 



are considered [^, |21|, |2^. We shall restrict our study of interactions here, however, to 
the case of the U(l) electromagnetic gauge. 

The general outline of this work is as follows. In Section II we present a short 
summary of the basics of the formalism we use. We then establish our basic requirements 
of the equation, and the logic in deriving it. In Section III we define the basic structure 
of the equations of motion, the form of a continuity equation, currents, and probability 
density. The basic requirements force us to seek nilpotent matrices as building blocks for 
the equation, so we investigate all possible nilpotent 4x4 matrices in Appendix ^ (and 
state some useful facts about them), and we state the Lorentz covariant appropriate 
nilpotent forms. We observe that it is not possible to construct a single Schrodinger- 
Stueckelberg type equation by iteration. We manage however, in Section IV, to construct 
two equations of motion, for transverse and longitudinal modes. After doing that, we 
still have four possible forms for both versions, and we show the relations among them. 
Two of them vary by the relative signs of their chiral components, and the other two are 
complementary in a sense that is described later on. In Section V, we investigate the 
form of the resulting currents and establish the validity of the form of the probability 
density for both versions. At this stage we observe the connection between the transverse 
and longitudinal equations, and the spin and convection currents, respectively, of the 
Gordon decomposition of the Dirac current. In Section VI we define the tensor product 
Hilbert space over the longitudinal and transverse modes, and show that our formulation 
is fully compatible with the one of Horwitz and Arshansky [Q. Next, in Section VII, we 
introduce a minimal coupling for the gauge field, producing electromagnetic interaction. 
We discuss the second order form of the generator of motion, and show it coincides with 
the one achieved by Horwitz and Arshansky EM. It implies the correct gyromagnetic 
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ratio, and is fully Hermitian under the positive definite, invariant scalar product for 
the quantum mechanical Hilbert space. In Section VIII, we show that we recover the 
original Dirac equation and its solutions in the limit where is constrained to be 
proportional to on mass shell. The solutions of the equations of motion are discussed 
in Appendix 0. We explore the discrete symmetries of the theory, in Appendix 0, and 
define the generalized parity, charge conjugation, and r reversal transformations. The 
states that are transformed one into the other under these transformations are pure left 
or right handed. Finally, in Appendix we make explicit the representation of the 
Dirac matrices we used for calculations, which lends itself to the problem, the chiral 
representation. 

We use the notation 

Xl-t ^ ( /vjl /y2 /vi3\ 

and we shall always use the symbol as the time component, and reserve for the 



time-like vector. We use an opposite metric relative to the Bjorken-Drell convention 

fl'oo = — 1 = —Qii = —922 = —933 
Furthermore, we shall use the notation 

where is a four-vector. We also use the uppercase P^j to denote an operator, and 
lowercase for eigenvalues. 

II. r-formalism 

We refer to the invariant parameter as r, the invariant universal world time, which 
describes the evolution of an event moving through space-time. The notion of the need 
of some invariant evolution parameter to replace the covariant time is not new, and we 



give some of the arguments as presented in [25 



Non-relativistic quantum mechanics uses the Newtonian universal time to describe 
a state in terms of square integrable functions over three-space at a specific time, which 
evolve according to Schrodinger's equation. But the Hilbert spaces associated with 
different times are distinct; we cannot superpose wave functions at different times. This 
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situation is inconsistent with special relativity. Viewed from some other frame the wave 
function becomes a function on different times, thus losing its interpretation as a state. 

The Klein-Gordon and Dirac equations have been able to resolve the problems of 
covariance; they have a manifestly covariant form. But the problem of constructing 
localized states still remains. The inconsistency of the solutions of the Klein-Gordon 
equation as an amplitude for a local probability density was shown by Newton and 
Wigner [||]. They showed that the distribution corresponding to a localized particle is 
an eigenfunction of the operator 



and the wave function corresponding to a localized particle is spread out by the order of a 
Compton wavelength. They reached similar conclusions concerning the Dirac equation. 



support, or localized in some other sense, does not maintain its localizability, and evolves 
out of the light cone, i.e. acausally. 

Quantum field theories make the transformation laws of special relativity and quan- 
tum mechanics consistent by assigning the spatial variables to the same parametric role 
as the time. The dynamical variables are operator valued fields which are functions on 
this parametric space-time. But the one particle sector wave functions of such theories 
describing the transition amplitudes between the vacuum and one particle states, suffers 
from the same difficulties as mentioned above. Predictions of phenomena concerning 
local properties in space-time are very difficult to formulate and interpret (e.g. interfer- 
ence phenomena), although spectral properties of non-local observables (energy. Lamb 
shift, anomalous magnetic moment), can be computed and are in excellent agreement 
with experiment. 

The basic difficulty of developing a consistent theory which incorporates the ideas 
of special relativity and quantum mechanics is related to questions concerning the re- 
lationship between time and locality. On one hand, space and time transform with the 
Lorentz group (and thus this relativistic time has a geometrical interpretation), and on 
the other hand, in a specific frame, it has been considered as a measure of evolution, 
of change. A way to resolve this ambiguity is to define the state of a system in terms 
of a distribution of events in space and time, while their evolution is parameterized by 
the time indicated on an ideal clock which is associated with every inertial frame (see 




(2) 



Hegerfeldt [Q, has shown that a distribution at a specific time, defined to have compact 
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p7| for further discussion). We call this parameter the universal time r, and it can be 
identified with Newton's time. 

Stueckelberg |Q , Horwitz and Piron , and others [0, ^ |10[ developed an underlying 



formalism incorporating the invariant parameter r which enabled them to construct a 
consistent manifestly covariant relativistic classical and quantum theory. The Newton- 
Wigner as well as the Landau-Peierls |^ problems have been understood in this 
framework. Two body problems, both for bound states and scattering have been treated, 
and the Zeeman p9| effect and selection rules [0 for radiation worked out. The spinless 



theory, in interaction with radiation (U(l) gauge field) has been second quantized ||3l[] . 
We state briefly the main principles of the formalism. 

The equations of motion of the classical theory, may be derived from the Hamilton 
principle p[ 

S j{v,dq^-K{p^,q^)dT)=Q (3) 

This principle is equivalent to the canonical equations 

dpfj, OK 



dr ^q^^ 
dqi" dK 



(4) 



dr dp^ 

describing the motion of an "event" along its world line (trajectory). For example, for 
a free event one takes 

° 2M 2M ^ ^ 

where M is a given property of the event, and sets the scale between r and the quantities 
of motion. We then have 

p dt E 

d^ ~ M ' d^ ~ M ^ ^ 

The proper time interval for the motion of a free event is defined by 

ds^ = de - rfx^ (7) 

and satisfies ^ 

ds' = ^,dr- (8) 

where m? = i?^ — is a dynamical variable to be determined by initial conditions 
and dynamics of the system. If initial conditions are chosen so that = M^, the 
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"on-shell" condition, then the proper time interval and the universal world time r, 
coincide. It should be noted that the theory is not constrained to time-like motion, so 
tachyonic propagation of events is possible. However, this does not imply explicitly the 
existence of tachyonic particles in laboratory measurements since one generally observes 
asymptotic states as the initial and final outcome of collision experiments. Due to 
the asymptotic conservation of the generator of free propagation for example, in 

potential scattering, states that are initially in time-like motion are time-like in the 
asymptotic final state as well. The structure of the theory, however, does not exclude 
tachyons a priori, classically or quantum mechanically. 

In the quantum domain, the states of the system for a given r are described in the 
Hilbert space L'^{R^,d^xdt), the space of square integrable functions of four variables 
given for the spinless 

(V") x) = 1 4'*xd^^- We do not define the scalar product 
for spin | particles at this stage; instead we shall show that we obtain the scalar prod- 
uct suggested by Horwitz and Arshansky p3|, and Arensburg and Horwitz |^, from 



considerations of the equation of continuity. 

The observables of spacetime coordinates and momenta satisfy the commutation 
relations 

^[P^Q1=^7'^n (9) 
The evolution of a state vector is described by the Schrodinger type equation 

idriJ = K%Ij (10) 

where for a free particle K = resulting in the Schrodinger-Stueckelberg equation 

(ID- 

Horwitz and Arshansky have suggested a second order equation for particles with 
spin I . They argue as follows. For a particle with spin the components of the wave 
function must transform as a representation of the Lorentz group. The norm must be 
invariant, so the representation must be unitary. But the Lorentz group is a non-compact 
group, therefore the unitary representations are infinite dimensional, containing all spins; 
such a ladder representation would, however, introduce problems with the application 
of the Pauli principle, for example, in the Sommerfeld model of a metal. If one were to 
use an induced representation based on the particle four-momentum as done by Wigner 
33| , the expectation value of x^, which by Eq. is replaced by i-^, would not be 
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covariant (the derivative acts on the unitary operator of the httle group defined by p^). 
Their solution consisted of introducing a representation induced on the httle group of 
a unit time-like vector, which we denote here by t^, which commutes with and p^. 
They described the transformation properties of the wave function, and found the form 
of the positive definite, covariant norm to be 

N = [ d^X^rt{x){^-t)iJrt{x) (11) 



where iprt{x) is the Dirac spinor, and 7'^ the usual Dirac matrices. They constructed 
the Hermitian and anti-hermitian parts of the operator (7 ■ P) under the scalar product 
associated with the norm. These are, in Hermitian form (under the norm (pT])), 



Kl = l{{i-P) + h-t){^-P){i-t)) = -{P.t){^-t) (12) 
Kt = ]p\{^-P)-{^-t){^-P){^-t)) = -2i^\P-K){^-t) (13) 

where = S^^tj, , T,^'^ = and the subscripts T and L denote transverse and 

longitudinal parts relative to the time-like vector t^. Since 

Kl = {P- tf ■ K^ = P^ + {P- tf (14) 

for the equation of evolution, Eq. (0), one can write 

idri^ = ^{Kl - KDi^ (15) 

By introducing minimal coupling — > — eA^ they then obtained 



where 

= S''^ + KH"" - K'^t" (17) 

This equation reproduces the correct gyromagnetic ratio, and does not contain the non- 
hermitian spin term appearing in Dirac's electromagnetically coupled, second order equa- 
tion (in the special frame = (1,0,0,0), one easily sees that S^-' is canceled so that 
there is no direct coupling of the electric field with spin in this special frame). Note 
that T,^^ F^j_y, appearing in the Dirac second order equation, contains io" ■ E as well as 
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cr ■ H; the former is not Hermitian in the norm / ilj*{x)ilj{x)d^x, which is Dirac's scalar 
product. 

In a later work, Arensburg and Horwitz p2| extended the formalism to a first order 



equation for spin | . Since the Kx part is responsible for the production of the correct 
gyromagnetic ratio, they postulated a first order equation of the form 

idrip = 2{P ■ KT){-f t)'ilj (18) 



(Their Kt is equivalent to that of Ref. |^ up to a factor — ^7^). Furthermore, they 



found its solutions and associated current, and showed that the current, although ex- 
hibiting a space-like nature at each point on the orbit (defined by t^), integrated over 
all possible in the forward light cone (completing the natural scalar product of an 
induced representation), reduces to a time-like current vector. 

Eq. (|18D does not, however, conform to the Schrodinger-Stueckelberg equation (|l|) 
by iteration; it leads to a second order equation in r. This is our original motivation in 
trying to obtain a new kind of first order equation for spin ^ particles in this framework. 



III. General features of the equations 
3.1 Basic structure 

In our attempt to find a manifestly covariant equation for an event with spin we require 
some general features of the equation concerning the relation between first and second 
order equations, namely the Dirac equation and the Klein-Gordon equation. We proceed 
from the point of view that a first order equation is an additional condition on the second 
order one, while at the same time introducing the notion of spin. In the process we 
narrow down the number of available options and maintain only the suitable ones. 

The desired equation should contain first order derivatives only, giving equal footing 
to the treatment of space and time. Another requirement is that the equation be Lorentz 
covariant, i.e., we wish the equation to be form covariant in respect to the choice of 
inertial frames. We assume that the space-time derivatives are coupled to an object 
constructed from 7^ matrices and perhaps some other four-vectors, and work with Dirac 
spinors, since we want the theory to be as close as possible to the standard theory of 
Dirac. For the actual Lorentz transformations of the wave functions and equations we 
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use the well known form |2^, 0] (this form is also applicable to the more general case 



treated in [|TT|, |^3|]; we use a'^'^ in place of Ti'^'^ henceforth for notational simplicity) 

^(A) = e-i^'^'^M^' (19) 

where 

^^^^ = ^[7^71 (20) 

and ujfj^u are the antisymmetric transformation parameters. 

We also require the equation to conform to a form of the Schrodinger-Stueckelberg 
equation by iteration, thus ensuring the free solutions of the spin particle to be also 
solutions of that form of Schrodinger-Stueckelberg equation, in parallel to the relation 
between Dirac's equation and the Klein-Gordon equation. This way each component 
of the wave function satisfies the free form of the Schrodinger-Stueckelberg equation 
separately. From another point of view we may regard the desired equation as an 
additional condition on the solutions of the Schrodinger-Stueckelberg equation, as the 
Dirac equation is to those of the Klein-Gordon equation. 

Now we postulate the most general form of a first order spin equation 

L{P)iIj = s.N^idrip + s+MN+ij (21) 

where L{P) is a linear function of first order space-time derivatives, namely a function 
of P^, and s+ and s_ are sign variables to be determined. and are unknown 
matrices at this stage. The second term on the R.H.S. of Eq. (|2T|) has no derivatives. To 
compensate for the dimensional deficiency we introduce the scale factor M, which appears 
in the Schrodinger-Stueckelberg equation. This term is necessary for maintaining the 
first order derivative in respect to r, after the iteration is done. 

To get the Schrodinger-Stueckelberg equation we multiply by L{P) on the left 

L'^{P)i(j = s-L{P)N^idri) + s+ML{P)N+i) (22) 

Now, let us define the commutators and anticommutators for the operators in Eq. (pl|) 

Ci^ ^ {^(P),iV_} 
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Ci^ ^ [L{P),N_] 

Cnn = [iV+,iV-] (23) 
We substitute Eq. (|21|) into Eq. (|2|) to obtain 

L\P)iIj = Nldli) - s+s.MCnnidr^ + s^QAdrt^ + s+MCi+i) - M'^Nli) (24) 

We do not want second order derivatives in respect to r, so we must have 

A^! = (25) 

Since we postulated that only first order space-time derivatives appear in L{P), and 
since Eq. ( p^ is supposed to coincide with the Schrodinger-Stueckelberg equation (|1|), 
we require in addition to Eq. ( pSj) also that 

= = (26) 

This gets rid of first order space-time derivatives in Eq. (p^, and leaves only second 
order ones in L^{P). We still have an unwanted term, M'^N^ip, and we require it to 
be zero (we could alternatively absorb it in the phase of the wave function). Thus, we 
obtain 

[-s+s.MCnn]-'L^iP)^(c) = idr^^c) (27) 
which is of the form of the desired equation (|1]), if we can find a solution for the conditions 

L\P) = ; = ±2 (28) 

The subscript (c) in Eq. (|27| ) represents an equation derived using anticommutators. 

We may use commutators instead of anticommutators, and repeating the previous 
procedure we reach the same conclusions concerning iV+, and Cnn- However, this 
time there is the requirement that 

Ci+ = = (29) 
The equation analogous to Eq. (^) is 

[s+s_MC'„„]-1l2(P)^(,) = zS.V'Cc) (30) 

The subscript (c) represents an equation derived using commutators. As we see later on, 
we shall need both forms. In addition, we require that the equation have a continuity 
equation, conserved positive definite probability, and currents. 
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3.2 Continuity equation and currents 

We are interested in achieving a continuity equation for Eq. (^) of the form 

9„r = ; (« = 0,l,2,3,r) (31) 

in order to obtain the usual interpretation, where j"^ = p should be the probability 
density, conserved through Eq. (|31|). The simple way to do this is clearly shown in all 



textbooks [0, ^ for the Dirac equation, and we follow this general method. We 
observe that since is nilpotent, the probability density cannot be of the usual Dirac 
form ilj^^ip, but contains a matrix between ip'^ and ip (a nilpotent does not have an inverse 
so there is no way to get a pure idrip term in Eq. (|2l|) ). Refs. |3^ show that the 
scalar product for the second order equation obtained from Eq. (|TT| ) is 



\4' Diraci "^Dirac) 

/ d^p (p Dirac (7 ■ t) l/j Dirac (32) 

where is the same time-like four-vector with norm = -1 mentioned earlier. The 
time-like nature of is crucial for the scalar product to be positive definite. [Of course, 
0, depend on f^, so that (|32D is a covariant structure on a bundle (i.e., as an induced 
representation).] The probability density is just the special case of the integrand of the 
scalar product when (p = tp. To find the continuity equation one multiplies the equation 
from the left by ip'^'j^ to get 

ijL{P)i) = s_tpN_idri^ + s+MtpN+tp (33) 

then multiply the conjugated equation from the right by 'y^ip to get 

^p^L\- P)-fy = -S-ipH Or Nl-f^^p + s+M^p^Nl-f^ij (34) 

and subtract Eq. ( p^ ) from Eq. (|33D. In order to achieve the form (|31| ) we must have 

7°(iV+)t70 = iv+ 

/(iV_)y = (35) 
7°(Lt(-P))^o ^ 

so that we get 

V'(L(P) + L{p))i^ = s_ij{N4i dr +1 dr])tp (36) 
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and finally 

d^iipb^ifj) = s-driipN.-^) (37) 

where it is assumed that the form of L{P) is b^idf^, and is some yet unspecified object 
with a Lorentz index which couples to P^. We must therefore now search for nilpotents 
for which {il)N_il)) is positive. A general analysis is carried in Appendix 0. 

3.3 Acceptable nilpotents 

One finds (in Appendix |^) that the candidates for nilpotents which pertain to a positive 
definite probability density are: 

(7-/W) : (7-t)(l±7') 

(7-/W) : {-i-t)±ats (38) 
(7 ■ t) ± Z7^ : (7 ■ t) ± ats 

where we denote 

ats = a^'^t^Sy (39) 

and Z, s, and t are light-like, space-like, and time-like vectors respectively. 

Nilpotents of this type come in non-equivalent pairs, in the sense that no unitary 
transformation connects the members of the pair (see Appendix 0). Furthermore, pairs 
in the columns of Eqs. (^) may be equivalent since we are able to transform from one 
to another through a unitary transformation. In any case we look at all four possibilities 
in the quest for L{P). Each nilpotent pair can be seen to be formed of two parts, a 
projection operator Pr±, and (7 ■ t), such that (7 ■ t) ■ Pr_ = Pr+ ■ (7 ■ t). In particular, 

(7-/W) = (7-t)(l±(7-t)(7-s)) 

(7-t)±z75 = (7.t)(l±z(7.t)75) 

(40) 

(7 ■ t) ± (7 ■ t)75 = (7.t)(l±75) 

(7-t)±ats = (7-t)(l±z(7-s)) 

The last of Eqs. (|^) is due to the fact that (t ■ s) = and ats = ^[(7 ■ t), (7 ■ s)]. 
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3.4 The four- momentum part : L{P) 

We wish L{P) to be linear in first derivatives, have zero anticommutators with N± (see 
Eq. (p6D ), produce the four-current part of the continuity equation, iterate to ±P^, 
and be Lorentz covariant. We check the four available forms of nilpotents in Eqs. (|38D , 
stemming from the specific reference frame in which they were found. To comply with 
Eq. (p6D we search for all F's (see Appendix ^ which anticommute with N±, then 
we couple in all possible ways, and check the anticommutators, and other criteria. 
The same procedure must be repeated for commutators of Eq. (^9]). Going through all 
options of Eq. ( ^OD is a tedious process, and we shall demonstrate only an example. Let 
us partly analyze the nilpotent form of N± = (7 ■ l^"^^). As the original non-covariant 
form we take 7*^ ± 7-*^. The F's that anticommute with 7" = F2 and 7"*^ = —iT^ are : 

F4,F5,F6,Fi6 ^ 27', ^7', 7V, 7' (41) 

Since it is required that L'^{P) must be a multiple of the identity matrix, if we wish to 
put in any combination of the matrices appearing in Eq. (^), we must divide them into 
two sets, for which the matrices in a set anticommute among themselves, and commute 
with all other matrices of the other set, and consider each set separately. Otherwise, 
after the iteration of L{P) we would find non-diagonal parts due to the anticommutator 
of two commuting matrices. The two sets are (7*^7^) and (27^,27^,7^). For example, the 
combination of 7*^7^ can be coupled to in a covariant way as 

7°7iPot;i ± 7°7^PiWo ^ (^Pv (42) 

where P^ = (Pq, Pi, 0, 0), and v = {vq, vi, 0, 0) is some vector, in this frame (since there 
are two matrices we form a tensor). If we want to conform to Eq. ( ^61) by taking (7 ■ Z*^^)) 
as A^_|_, and we do not want P^ to be only light-like, then v must be light-like (since 
= {aPv, 7 ■ = implies P oc or f oc Z^^*"^). But then the vector /^^^ must be 
the same in both N^,N_ (either l^^^ or /*•"■*), which makes = 0. This contradicts 
our basic assumptions, so this sub-form must be discarded. 

We have required in Section III that the iteration of the first order equation should 
give us L'^{P) = ±P^. Checking all available options of the nilpotents of Eq. (^0]), indi- 
cates that we are unable to get such a form by iteration (at least not under our initial 
assumptions concerning the iteration process). [In Section VIII we shall show how to 
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produce by a somewhat different iteration procedure (obtained in a rigid limiting 
case); in that procedure we lose the Hermiticity of the spin term in the second order 
electromagnetically coupled equation, getting an expression as in the Dirac equation. 
Furthermore, we lose the structure of the continuity equation.] This applies to commu- 
tators and anticommutators alike. Nevertheless we proceed, considering the transverse 
and longitudinal options separately, combining them in Section VI. We start by choosing 
from the anticommutator forms. From them we choose 

7V± = (^.t)(l±^5) (43) 

and 

L±iP) = a^'^P.U (44) 

This is because in doing so, we deal with only one additional four-vector, t^, while 
retaining the ability to achieve the gyromagnetic ratio, and exhibiting some features 
concerning chirahty due to the use of the projection operators |(1 ± 7^). It also seems 
that this is the simplest choice which has as many benefits as possible in this situation. 
As we show later on, this choice gives rise to the transverse Schrodinger-Stueckelberg 
equation. Not excluding the existence of a longitudinal equation also, we follow our 
current choice for N±, and take the simplest choice for the longitudinal equation, namely 

L\\(P) ^ i(P ■ t) (45) 

IV. The equations of motion 

We now discuss the consequences of the above choice for the forms of N± and L(P)'s. 
This discussion concerns what we call the transverse and longitudinal Schrodinger- 
Stueckelberg equations. 

4.1 Transverse equations of motion 

We have 

Cnn = {(7 • ^)(1 - 7'), (7 • + 7')} = -4^^ (46) 

taking = — 1, a unit time-like four- vector, we have to arrange an additional factor of 
2 to get the desired form of We place this factor with iV_|_. [In Section VIII we show 
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how to get the solutions of the Dirac equation in the hmit where and coincide. For 
that to happen we have to take the factor 2 as indicated.] The anticommutators satisfy 
of course 

{(7-t)(l±7'),(aPt)} = (47) 
which follows from the important relation 

{(7-t),(aPt)} = (48) 

Iterating Lj_{P) we obtain 

{aPtf = Ph^ - (P ■ tf (49) 
We define the transverse and longitudinal momentum, relative to t^, as 

Pl^. = P^ + {P-t)t^ (50) 
P\\, = -{P-t)t, (51) 



and we therefore have 



The equation of motion is 



[aPtf = t^Pl = -p2 



- hc{a^'id^K)ij^ = hs_{j ■ t)(l - ^')tdriJ± + s+^i7 ■t){l + 7^)^^ (52) 

We display h and c in this principal equation; elsewhere we take h = c = 1. It is 
understood that x'^ = ct. Notice that in our metric P^ = —i-^- 
Iterating, we obtain 

= idriJ± 

s+s_2M 

This forces the signs (with our choice of phase) to be either both positive or both 
negative, so that = 1, and we get the transverse Schrodinger-Stueckelberg equation 

— —i'± = idr'ip± (53) 
2M ^ ^ 

which we assume to hold. Using the notation of Horwitz and Arshansky , 

L^(P) = —K^, which generates the evolution of transverse momentum only. Note that 

iIj± contains only the transverse space-time coordinates. 
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We have not excluded the possibihty of a longitudinal Schrodinger-Stueckelberg equa- 
tion; from another, independent, first order equation, which we shall discuss later, we 
shall also have 

^^11 = ^^r^p\\ (54) 

where ipw contains only longitudinal space-time coordinates. The subscripts _L and || are 
introduced for making a distinction between the wave functions of the transverse and 
longitudinal equations. These indeed constitute kinematically independent degrees of 
freedom. 

Since Eq. ([5^) contains projectors we can decompose it into two coupled equations 
by multiplying it from the left with the same projectors. Denoting the wave function as 
composed of two chiral spinors, (two-component spinors in the chiral representation), 

V'L = ^(l-7')^ V'R = ^(1 + 7')V' (55) 

we obtain 

{^■t){aPt)^^L = s+M^^R (56) 

(7-t)(aPt)V'±i? = 2s.idriJ±L (57) 

This is an explicit chiral decomposition, not symmetric for left and right handed spinors, 
a matter discussed later on. Such a form for equations of the first order in r-formalism 



has been proposed by Davidon [Q, in a somewhat ad hoc manner. He could not over- 
come the problem of finding a positive definite probability density, and he obtained 
a non-hermitian spin term after coupling the electromagnetic field. The equations 
(|56D,(|57D are completely equivalent to the one equation (0); it is just a rewriting. 

Concerning the sign variables s+ and s_, it seems that we have two choices. Denoting 
by 'ipi_+ and ip^_ the solutions of the two sign versions of Eq. (|52|) , by adding and 
subtracting the two versions, we get exactly the equations (^6D,(p7D, with the following 
relation between V'x+^V'x- and 4'±l,4'±r 

+ ij±- = ij^R (58) 
V^±+ - = (59) 
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As can be seen from Eq. (pTj), the sign s_ (equal to s+ in our treatment), sets the sign 
of the so called Qrr component of the metric. 

After we have discussed the two options of sign assignment, we cannot ignore the 
two possibilities of nilpotent assignment. Instead of deducing Eq. (^2|), we could equally 
well introduce a complementary equation with the places of the nilpotents interchanged 

- a'^'td^t.ips. = 5_(7 ■ t){l + 7')^5,<^x + s+y (7 ■ - 7')<^± (60) 

where we left the signs in place. Eq. ( |60D has all the features of Eq. (|52D, except that 
the role of the left and right handed spinors is interchanged. Writing the two relevant 
coupled equations analogous to Eqs. (p^),(p7D we get 

{l-t){aPt)ip^R = s+Mif^L (61) 
(7-t)(aPt)(/.XL = 2s_idr^^n (62) 

Now, if we abide by Eq. (|52D, and describe the theory as portraying the evolution of 
a left-handed spinor, ip_LL, while ip±ji is just a sort of auxiliary field, we give up right- 
handed events. This cannot be done without some justification. In Section we see that 
such an equation gives rise to a continuity equation, (through a simple procedure similar 
to the one performed for the Dirac equation), where the probability density is composed 
of left-handed spinors only. 

The interpretation of tp±[i as an auxiliary field, justified by Eqs. (|57D ,(|56D where only 
'ip±L is seen to evolve according to r, the existence of the other equally justified equation 
(|60D , the inability to prefer Eq. (|5^) over Eq. (|60|) , and the requirement for a reasonable 
probability density, lead us to a unification of the two. 

We define 

0± = i^±L + <^±R X± = i'^R + '^±L (63) 

These are now four-spinors, where x± is an auxiliary field, and is the main field, 
evolving with r. The two sets of coupled equations (|56D,(|57|) and (|61]) , (|62D , become 

(7-t)(aPt)0^ = s+Mx± (64) 
{^■t)iaPt)x± = 2s^tdr(l)± (65) 

by just adding them up. We denote them as the "extended" form. We therefore have 
achieved a satisfactory probability density of four-spinors, which is sufficient for the two 
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irreducible representations of SL{2,C) (see |2^,[^). Of course by iteration we get the 
transverse Schrodinger-Stueckelberg equation for (f)± and for x±- 



4.2 Longitudinal equations of motion 

In dealing with a first order equation that agrees with the longitudinal Schrodinger- 
Stueckelberg equation, Eq. (|5^), we check it for the nilpotents and momentum part we 
choose. Since the nilpotents are the same as for the transverse version all conclusions 
from Subsection apply here too. Iterating L\\{P) we get = —{P-ty, and the relation 
between s+ to s_ is maintained as in the transverse case. The equation of motion is 

t{-td,t'^)^\\ = s_(7 ■ t){l - l')tdr^\\ + s+^{^ ■ + 7')V'|| (66) 
Using projectors of Eq. (^) we can decompose it in two 

i{-i-t){P-t)iP\\L = s+M^\\R (67) 

■ t){P ■ t)ij\\R = 2s^tdr^\\L (68) 

The discussion concerning both options in using the sign variables applies here as well, 
and we do not repeat it. As for the transverse case, Eqs. (|6^),(|68D can be used to 
show that the left and right handed spinors pertain to the longitudinal Schrodinger- 
Stueckelberg equation. 

Exchanging the place of nilpotents in the equation to get the complementary equa- 
tions of motion for the right handed part as the main field, while maintaining the inter- 
pretation of the auxiliary field, and adding the two types of equations as in the transverse 
case, leads to the coupled equations 

z(7-t)(P-t)0|| = s+Mxw (69) 
z(7-t)(P-t)X|| = 2s^zdr<Pii (70) 

with the same structure of and x\\ as in the transverse case. 

V. Probability density and currents 

The continuity equations for the transverse and longitudinal versions are obtained by 
using the procedure for obtaining the primary form of a continuity equation in Subsection 
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. We show that the currents form consistent continuity equations, and are related to the 
Gordon decomposition of the Dirac current. 

5.1 Currents of the transverse equation 

Using the procedure in Subsection on Eq. (|52|), we obtain 

- d^i^P^a^^'^Uij^) = 2s.dr{Ml ■ - 7')^±) (71) 

We observe that the R.H.S. of Eq. ( [7TD is a probabihty density, exactly in the sense de- 
fined by Horwitz and Arshansky but achieved without exphcit group theoretical ar- 
guments. It is obviously positive definite for wave functions containing left handed com- 
ponents, which can be seen by transforming to a reference frame where = (1, 0, 0, 0). 
We denote this probability density (in the chiral representation of 7 matrices) as 

P±L ^ {Ml • ^)^(1 - 7')^±) = ±L{^%)i^±L (72) 

where 

a^^(l,-a') ; a^^{l,a') (73) 

and cr* are the Pauli matrices. The appearance of p±L only in Eq. ( |7ID is the chief 
reason for introducing both Eqs. (p^) , (pDj) . From Eq. (^) we obtain in the same way 
p±R, and we can add them together, and obtain = p^^ + P±r, thus having the 
a density containing both SL{2,C) representations as in [^, Since p±R is also 

positive definite, so is p±, as required. 

We define the four-currents in Eq. ([7T| ) as 

jlL = M^'"'tu)i^± (74) 

so we actually have 

- ^Mi±L = '^S_drp±L (75) 

which is the five dimensional conservation theorem required. Let us consider 

i^±L{crPt)^^R = 2s_i:^Lh-t)idri^^L (76) 
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and the conjugate of Eq. (|7|) 

Using the equations of motion, Eqs. (^6l),(^7|), we find 



(78) 



dri^LL = -2ir(7-i^)(^^^)^±R 

Using Eq. (|56|), and the sum of Eqs. ([76|), (|77D we get 

drP±L = -i^^LL{{c^Pt){l-t){aPt) + {^-t){aPt){aPt))^^L 

' Ml-t){-Pl + Pl)ij^L (79) 



2M 

To simphfy the expression for p±L we state some useful relations 



P±^ = ~id±^ = -idf, + {-idy ■ f 
In the reference frame where = (1, 0, 0, 0) the transverse derivative is seen to be 



Pi = -d±,d^^ (80) 



So, 



d±o = 
d±i = di 



drP±L = ^d^,{^^L{l ■ t) ^xl) (82) 



Using Eqs. ([Tq) we find that the r derivative of the auxiliary field also obeys a continuity 
equation, exactly as the main field with ip^^ exchanged with ip^R 

i - ^ 

drP±Lau. = ^5±/.(V^±ij(7 ■ t) dt_ i^±R) (83) 
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where we therefore understand p±Laux a "density" associated with the auxihary field 
of il'±L, namely with ip±R. This form of current resembles the Klein- Gordon current, and 
when taken in a frame where = (1, 0, 0, 0) we recover the space part of the current. 

To make a connection with the form of currents exhibited by the Dirac equation, we 
re-express Eq. ( fTll) as 

-UjlL = 2s.t,drplL (84) 

where 

jIl = d.i^P^a'^^i;^) (85) 
pIl = ^^lY^^l (86) 

Now we integrate over a space-like hypersurface enclosing a four-dimensional volume. 



where is normal to the surface (see, for example, ||3^). Therefore, we imagine a 
space-like hypersurface perpendicular to and another one like it some distance above, 
meeting it only at infinity; it then follows from Eq. (p^ ) that 

- / c/S KjIl = 2s-dr I dL t^p'iL (87) 
where denotes the closed space-like hypersurface. We transform to the four- divergence 

- j dv BJjIl = 2s-dr I dv d^p^L (88) 

JVt JVt 

But since 

- = 9,9^(^xa'^'^^^) = (89) 

on account of the antisymmetry of a^^ and the symmetry of d^d^^ we get 

/ dvdrd^p^^L = ^ (90) 

JVt 

Now, assuming this is true for the integrand also, i.e., for any thin slice V^, we obtain 

drd^ii^^LYij^L) = (91) 

and hence 

dM±Lri^±L) = f{xn (92) 

The quantity d^i^jj^iY'^LL) is independent of r, and if it is not zero, it is some r- 
independent function f{x^). 
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Comparing to the Gordon decomposition of the Dirac currents one finds a similar 
structure when applying the continuity equation. 

We now consider the combination of the above methods and equations concerning 
Eq. (^2|), with their application to the nilpotent interchanged equation, Eq. (^). The 
same line of thought, when applied to Eq. (|60D , yields the same results, only ip^ji replaces 
i^LL, replaces '?/'±_r, drp^R replaces drP^i- Ultimately we get the continuity equation 
for the main field 

dr{4>L{l ■ m±) = ^d^.iMl ■ t) ^ 0±) (93) 

Furthermore, one obtains for the total main field, a r independent continuity equation 
resembling the one in Dirac's theory 

5r5^(0±7^0±) = (94) 

If 9^(0±7^0±) = 0, (i.e., in case the sum of divergences for and ip^^R cancel), we 
have the form of the Dirac continuity equation, but for the transverse field (recall that 
the wave function 0^ itself still has dependence on t^). 

We have a physical behavior of the particles described by our equations, which re- 
sembles the physical behavior of Dirac's particles, without the use of Dirac's equation, 
but using instead the transverse equation. It is interesting to observe that the Dirac 
type currents were obtained from the probability density of our theory. We may view 
this as though our theory is in some sense a fundamental underlying structure for what 
is essentially Dirac's theory. 

We remark, moreover, that we have obtained in our method, a scalar product which 
agrees with that of 123], i.e.. 



d^x 0x1 (7 ■ t)(p±2 (95) 



5.2 Currents of the longitudinal equation 

Using the procedure in Subsection on Eq. (pB|), we obtain 



^11 a^r^ii =2s_za.(^||(7-t)^(l-7')V^I|) (96) 

In the R.H.S. we identify a positive definite quantity. Unfortunately, we cannot under- 
stand Eq. (pB) as a continuity equation. However, we can overcome the problem by 
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observing similar features of the equations as in the transverse case. We state some 
features of the longitudinal momentum 



\\ 

(97) 



and in the special frame where = (1, 0, 0, 0) we have 

<9||o = <9o 

d\\i = (98) 
For a treatment analogous to that of the transverse equation we find 

5rP||L = ^%(^||L(7-t)#^i|L) (99) 

Using the space-like hypersurface integration while defining 

J\\L = V^ii^^V'ii 

tjjl^ = 2s.t^idrP^^^ (100) 



we get 

/ dv2s_drd^{lP\\Lr^\\L)= - f dv d.iijp&'iJw) (101) 
JVt JVt 

Eq. ( p.01| ) is similar to Eq. (|88| ) in the sense that the probability density gives rise 
to a form of a Dirac current. The transverse current of Eq. ( pS]) has the form of the 
spin current of the Gordon decomposition of the Dirac current, while the longitudinal 



current of Eq. ( 101 ) has the form of the convection part. It seems that the two parts of 
the Gordon decomposition manifest themselves in the two versions of the equations of 
motion. In case <9^(^||l7'^V'||l) is independent of r, one would have a result for 
dy{ip\\id^ip\\), similar to that of the transverse part. If we integrate Eq. ( |101| ) over r, one 
obtains a conservation law of the same type as well. 

Finally we state the equation for 0||, the four-spinor main field 

^rP|| = ^%(0||(7-t)#0||) (102) 
The solution of the equations of motion is given in Appendix |B|. 
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VI. Product Hilbert space 



The vector splits the evolution of the momentum into two parts, the transverse and the 
longitudinal. These two modes of the motion of an event, transverse and longitudinal, 
are complementary. Therefore, the overall Hilbert space is the tensor product of the 
transverse and longitudinal Hilbert spaces 

n = n±® Hw (103) 

Considering the free event, the probability density is 

P = P± • P\\ (104) 

By defining 

= {^■t){aPt) (105) 

L|| = i(^-t)(P-t) (106) 

a; = 2M (107) 

and absorbing the signs s+,s_ into cp or x, the transverse and longitudinal equations 
can be expressed in a similar manner 

L±(p± = ujx± (108) 

L±X± = idr(l>± (109) 

and 

= ujx\\ (110) 

^IIXII = idrcpw (111) 
In the product Hilbert space we define 

= (112) 

X = X±®X|| (113) 
and L±,L\\ operate on their respective factors. Now 

L±<f) = uj{x±^(f)\\) (114) 

L||0 = uj{(f)^^x\\) (115) 

L±x = iidr(t)±) (E) x\\) (116) 

L\\X = X±®(^5.0||) (117) 
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Denoting 



we obtain 



and 



L = L±® L\\ (118) 

L0 = uj\ (119) 
Lx = {idr(t>A.)®{idr(t>\{) (120) 

(121) 

= L^x = {(^idr(p±) {(jJidrCpw) (122) 



We are interested in getting a unified equation for the transverse and longitudinal 
modes, conforming to the full Schrodinger-Stueckelberg equation. Therefore we consider 
the combinations 

(L^ + L||)0 = a;(x±®0|| +0±®X||) (123) 
{L^ + L\\)x = {idr(l>±)^X\\+X±^{idr(f>\\) (124) 

Considering {L± + L^^Y we find 

(Lx + L||)V = u;idr(l) + 2u;'^x 

= ujidr(f) + 2L(f) (125) 

Since 

(L^ + L||)2 = Li + Lj + 2L (126) 

we obtain 

{Ll + L|)(/) = (127) 

which is the desired full Schrodinger-Stueckelberg equation. Repeating the same steps 
for the auxiliary field x, we find that it also obeys the full Schrodinger-Stueckelberg 
equation 

(Li + Lhx = ujtdrX (128) 
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VII. Interacting charged equations 

To introduce the electromagnetic coupling, we define the electromagnetic field appropri- 
ately on the manifold of the tensor product space i.e., 

a"{x,T) = a'^{x±,xii,T) (129) 

where a — 0,1, 2,3, t. We may therefore define 

= (^■t)(am) (130) 
Lf = i{l-m-t) (131) 

where 

Note that although L± and Ly act on the tensor product space as Lj^ (g) 1 and 1 (8) Ly, 
this is no longer true, in general, for and L^. Let us start by defining equations 
containing fields of the form a±n and ay^u 

(7 ■ t)K''(P^ - ea±^{x^, t))U]{Mx±, r) ® 1) = u;{x±{x±, r) ® 1) (133) 

■ t)[a^'''{P^-ea^^{x^,T))U]{x±{x±,T)^l) = [(i9^+eaxr(a;±, t))(/)x(xx, r)]®l (134) 

and 

i(7 ■ t)[(P^ - ea||^(a;||,r))t^](l ® 0||(a;||, r)) = cu(l ® X||(a;±,r)) (135) 

i(7-t)[(P^-ea||^(a;||,T))t^] (10X11 (^11,^)) = l®[{tdr + ea\\ r(x||, r))0||(x||, r)] (136) 
Now, since we can write 

53(al0x®a|^0||)(a;x,a;||) = J2<i^±miH)M^±)M^) (137) 
where n is some index, the limit of this sum can approximate any function a{x) i.e., 

'^a^{x±)a'\\{xii)4>±{x±)4>ii{x\\) = a{x)4>±{x±)4>\\{x\\) (138) 
Since we have gauge invariance 

a±n{x±) a±f,{x±) + ^d±f,A±{x±) (139) 
a\\^,{x\\) a||^(a;||) + ■^a||/,A||(a;||) (140) 
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we may generalize the gauge invariance to 

a^{x±,xi\,T) af,{x±,x\\,T) + ^d^A{xj_,x\\,T) (141) 

and 

ar(x±, x\\,t) — > ar(x±, x\\,t) drA(x±, x\\,t) (142) 

e 

These hold for the equations 

(7-i)((7m)^((/)^®0||) = u;(x±®0||) (143) 

{-ft){am)±{x±®X\\) = (i9^0±)®X|| + ea^(0±®X||) (144) 

2(7-t)(n-t)||(0x®0||) = o^(0±®X||) (145) 

i(7-t)(n-t)||(x±®X||) = X±® (^5r0||) + ea^(x±®0||) (146) 

since in each equation we can refer to the other variable (say || in the ± equation), as a 
parameter. The subscript on the matrix operators on the L.H.S. of the equations means 
that we operate with derivatives on the relevant factor space. 
We can therefore define the form of the basic equations as 

Ll<P± = u;x± (147) 

Ltx± = {idr + ear)(l>± (148) 

and 

Lf0|, = uxii (149) 

Lfxw = (^a, + ea,)0|| (150) 

Denoting as the tensor product operator obtained from and L^, as 

= Lf (151) 

we obtain 

L^c/) = co'x (152) 
L^X — {idT(f)±) {idr(j)\\) + e^Qr^cl) + earidrcj) (153) 
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As before, we obtain from the squared expression 

{L^ + L^f^ = uj{idr + ear)^ + 2uj^x 

= uj{idr + ear)(p + 2L^(l) (154) 

Therefore we find 

[(L^)' + (Lf )']0 = u{idr + ea,)0 (155) 

which is exactly the full electromagnetically coupled Schrodinger-Stueckelberg equation. 
Now, using the fact that 

{a^''n^U)iaf"'Uptx) = ^{[U^,Up][a^'''U, a^x] + {U^,Up}{af'''U, af^n^}) (156) 
the definition (this quantity was defined by Horwitz and Arshansky |^ as well) 

- 2ia^'' = [a^'^U, aPhx] = 2i{am'^ - {a^'Hxy - i(j''''U)t^) (157) 

and 

{a'^^'Vu, (tP^vx} = 2{v^g^'P - v^v^) (158) 



and the relation 
we get 
by defining 
Furthermore, we have 



[n^,n,] = ze^, (159) 



(L'lr = ni - ^arUu (160) 



= + (n ■ t)t^ (lei) 

{LfY = -(n-t)2 
ul = n2 + (n-t)2 



and it is easy to see that Eq. (|155|) is 



{zd^ + ea.)0 = ^[n' - l^rf,u]<P (162) 
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which is exactly the same equation found by Horwitz and Arshansky |2^. This shows 
that our formulation is fully consistent with theirs. 

The gyromagnetic ratio is taken as the ratio between the coefficients of the terms 11^ 
and |cr("^/^j,, giving the correct relation of relative size and sign, between the momentum, 
spin coupling, and mass terms (see |3^). Furthermore, we obtain, as in Horwitz and 
Arshansky a fully Hermitian spin term, as opposed to the Dirac case (and the 
additional field which they did not use). 



VIII. The Dirac limit 

We have seen that it is impossible under the assumptions of Section to find a single first 
order equation which iterates to the full Schrodinger-Stueckelberg equation. Therefore, 
we resorted to the combined product Hilbert space description in Section . However, by 
making a slight modification we can obtain the full Schrodinger-Stueckelberg equation 
by iteration from a single equation, but there is a price to pay; the gyromagnetic term 
in the electromagnetically coupled version is no longer Hermitian, like in Dirac's three 
dimensional case, and the structure of the continuity equation is lost. This modification 
is the necessary ingredient to see the connection between the theory developed so far 
and Dirac's theory. 

The modified approach that we use to employ a single equation is to break up the 
operator L{P) to two parts, which have different algebraic structure with respect to the 
nilpotents N±. This will enable us to obtain the full Schrodinger-Stueckelberg equation 
upon iteration. 

Let us assume that L{P) is composed explicitly of two parts, Li{P) and (P). 
Then, the general form of the first order equation is 

(Li(P) + L2{P))ip = s.N.idrip + s+MN+^jj (163) 

We now multiply Eq. ( |163|) from the left by (Li(P) — L2{P)). To be able to insert Eq. 
( |163| ) into the new equation we must have 

{Li(P),iV±} = 
[L2(P),Ar±] = (164) 
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Requiring that 

[Li(P),L2(P)] = (165) 

and that N± are nilpotents, we obtain 

(LliP) - LliP))^ = -s.s+MCnnidri^ (166) 

Performing an evaluation process similar to the one done in Section , we make the 
simplest choice which is 

U{P) = {aPt) ■ L2{P)^t{P-t) (167) 

and the nilpotents are N± = (7 ■ t)(l±7^). This is a synthesis of the use of commutators 
and anticommutators of Section , and of the transverse and longitudinal equations. The 
equation of motion is 

i-a^-'td^U + ^i-^^^t''))^ = s.{^ ■ t){l - ^')tdr^ + 5+^(7 ■t){l + 7^)V^ (168) 
Iterating we get the full Schrodinger-Stueckelberg equation (|l|). We obtain the currents 

- id,{^lja'%^lj) +ijd%ilj) = 2s_z9,(^i(7 ■ t)^L) (169) 
which can be written as 

-t^jt = 29,s_t^(^i7'^^i) (170) 

where 

j1 = (d^itfja'^^'tlj) + tfjidf"-^) (171) 

It is clearly seen that resembles the definition of the Gordon decomposition of currents 
of the Dirac equation. 

An equation with the roles of N± interchanged can built, resulting in a similar Gordon 
decomposition but this time only right-right terms appear in the current. When we 
discussed the transverse and longitudinal equations we could transform between the two 
versions opposite in the nilpotent assignment by a suitable transformation. Decomposing 
Eq. ( |168| ) by projection into two coupled equations we find 

-z(7-P)V'L = s+M^Pr 

-i{^-P)^R = 2s.idr^L (172) 
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At this stage the dependence of the equations disappears; we deal with this below. 
Taking these equations to mass shell, choosing s+ = s_ = — 1, and observing the 
eigenvalues of plane wave solutions (see Appendix |B|), we obtain 

-ti-fP)ijR = mijL (173) 
and for the equation with N± interchanged 

- 2(7 ■ P)>^R = -rmpL 

-2(7 ■ P)ipL = rmpR (174) 

To transform from Eq. (|173| ) to Eq. ( |174| ), we can multiply iphy —^{j ■ P) (working on 
shell) to exchange the roles of the left and right handed spinors, therefore the relation 
between and if is 

-yjR = ; Vl = (175) 



m m 



This implies that 



V^L = -#L ; = ii'R (176) 

and the relation between the main field and the auxiliary field x is 

(j) = tX (177) 

Using Eq. ( |176| ), Eqs. ( |173| ),( p!7^ ) are then exactly the Dirac equation in the chiral 
representation (holding for both the main and auxiliary fields) 

(7 ■ P)(f)R = -mcpL 

{7-P)^L = -m(t)R (178) 

This implies, of course, the Dirac current and continuity equations for the main and 
auxiliary fields (f) and x- 

When electromagnetic coupling is introduced in Eq. ( |168| ), assuming that a,- = 
(Coulomb- like gauge) and independent of r (restricting ourselves to the zero mode) 
for simplicity, we obtain for the second order equation 

^^.^ = ^[n^- (179) 
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which has the usual non-hermitian term, and conforms to the second order charged 
Dirac equation (when idr "H^? see [0). We therefore see that the decomposition 
into transverse and longitudinal modes was essential to achieve a Hermitian interaction 
as well as the Schrodinger-Stueckelberg form. 

All the discussion above can be approached from a different point of view. If we take 
the longitudinal equations of motion, and assume that the momentum is in the direction 
of the four- vector t^, such that 

P^ = at^ ; P^ = -a^ (180) 
then the longitudinal equations of motion become 

-z(7-P)V^I|L = s+M^\\R 

-z(7-P)^l|ij = 2s_z9.^||i (181) 

with a similar result for the nilpotent interchanged equation. Making the identification 
ip = ip\\, these are exactly Eqs. ( |172|) . One can now understand how disappeared 



previously in Eqs. (|172|) . Furthermore, looking at the basic structure of the solutions of 
the longitudinal equation for 0|| (see Appendix 0), Eq. ( p44| ) , and replacing by ^ in 
C|~(t), (on shell condition), we obtain the exact non-normalized solutions of the Dirac 
equation in the chiral representation 

where ^1^2 are two independent two-spinors; this is also apparent from Eq. ( |175| ) (actually, 
when taken in the energy representation of the Dirac matrices these solutions correspond 
to the positive energy, and C||''(^) solves the Dirac equation for the negative energy). 
The transformation of the extended helicity projection operator (see Appendix ^ is a 
bit dehcate since (crPt) and \ p± \ go to zero together. Observing Eq. (|168|) and its 
consequent Eqs. (|172|) , we find the extended helicity operator commutes with them. In 
Eqs. ( |172| ) does not appear, so we can take it to be = (1, 0, 0, 0), and the extended 
helicity operator becomes the usual helicity. We use the usual helicity operator to obtain 
four distinct solutions. 

On the other hand, for the transverse equation, taking in the direction of has 
different consequences. Observing Eq. (|5^), 



{a'"'P^P,)ij± = s_(7 ■ P)(l - 7')^5r^± + s+^(7 ■ ^')(1 + 7')^± (183) 
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then a^'^P^Py = 0, and p\ = 0. Considering the free solutions, there is no r evolution, 

= s+^(7-P)(l + 7')^± (184) 

and we can say that ip±R and (p±L are zero, there is no auxiliary field, and no transverse 
first order equations exist. All the evolution is described by the longitudinal equation. 

It seems that as departs from the direction of t^, we depart from the Dirac descrip- 
tion, and the event is described by two independent equations of motion, transverse and 
longitudinal. The event evolution in r becomes different for the two versions, differing 
in phase and spinor content. At the same time the auxiliary fields are spontaneously 
generated in the transverse case, and from the Dirac field in the longitudinal case. Fur- 
thermore, the auxiliary field, treated as a mathematical convenience, but seen to be 
strongly related to the main field, Eq. ( |177| ), may be a further indication of the amount 
of departure from the on shell Dirac theory. 

IX. Conclusions 

Accepting a second order spin equation of the form of the Schrodinger-Stueckelberg 
equation as the basic structure of the theory, and requiring that the free solutions of 
a first order equation be also the solutions of the second order one, lead to transverse 
and longitudinal equations. After a selection process we chose the most suitable forms 
for the ingredients of these equations. We introduced the induced representations on a 
time-like vector and an auxiliary field so that the theory be consistent. The evolution 
of the free event is governed by two complementary equations, the transverse and longi- 
tudinal. We manage to unify the complementary behavior of the event, transverse and 
longitudinal, with the use of a product Hilbert space. Then, after introducing the elec- 
tromagnetic coupling, we obtain the correct gyromagnetic ratio, and the second order 
charged equation describing the evolution of the main field, is fully Hermitian under the 
scalar product. Finally, we showed how the theory goes over to the Dirac theory in the 
limit in which is in the direction of on mass shell. In this limit the theory changes 
its structure. The first order transverse equation no longer exists; its auxiliary field and 
T evolution disappear, and the longitudinal equation takes on a new meaning, leading 
to the original Dirac equation. We also display, in Appendix 0, the solutions of the 
equations, which are characterized by the "extended helicity" and "chiral precedence" 
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for the transverse case, and by "extended helicity" and "extended parity" in the longi- 
tudinal case. Both cases have the same properties under the generalized parity, charge 
conjugation, and r reversal transformations (as shown in Appendix |C|). The objects 
which transform one into the other under these transformations are pure left or right 
handed spinors, thus exhibiting the chiral nature of the theory. 
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A Idempotents and Nilpotents 
A.l Idempotents 

Let us consider idempotents of the form (we carry out the analysis in a given Lorentz 
frame, and discuss later the corresponding invariant forms) 

(185) 

where are elements of the Dirac Clifford algebra (see Appendix 0). We state some 
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= 1 (186) 

rr' = r, (i87) 

Ti^ = Ti (2 = 1,..., 16) (188) 

where implies complex conjugate transpose, and for j ^ k, 

TjTk = ejkTi ejk e {1,-1, i,-i} (189) 
r^r, = iejk)-'Ti (190) 

There are fifteen pairs of Pj±. In fact they are equivalent under automorphism. We 
show this by defining a unitary transformation Tj such that 

TJ = ^{l + ^r,) ; T,-^ = -^(l-2r,)=r,t (191) 
Using Tj we can transform from one F to another. If Fj and Fj commute we have 

[F,-, Fi] = ^ [Tj, F,] = ^ TjT,Tj~' = F, (192) 
and because of Eqs. (|189|) , (|19CI|) we have (for i ^ j) 



[Vj, Fi] = ^ F,Fj - VjV, = {eij - e,r>)Ffc = 

^ = ^ = ±1 (193) 

and 

[F,,F,]=0 {T„Tj}^{] (194) 
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Fj is invariant under the transformation Tj; we find a triplet Fj, Tj, F^ which commute 
among themselves, where F^ = FjFj. These triplets form maximal commuting sets. 
On the other hand, if Fj and Tj do not commute then 

[r,-,r,]^o [r,-,r,]^o^(e,,-e-.^)rfe^o 

Sij ^ ej^ Sij = ±i (195) 
This means that for the non- commuting case we have 

[F,-,F,]=±2iFfc (196) 

and so 

Tjm-' = l(r, + i[r,-,r,] + F,T,r,) = .Frfe (i97) 

Furthermore we find that 

[F,,F,]^0 ^ {r,,F,} = (198) 

We can transform Fj into F^ using F^. This is a Pauli algebra structure for the triplet 
V v v 

'-it jt ^ k- 



[F,-,F,] = ±2zFfe 
[F,,Ffe]=±2iF,- 
[Ffe,F,-] = ±2iF, 



{r,,r,} = o 
{r,,Ffc} = o 

{F,,F,} = (199) 



From Pj-t and Pj± from a maximal commuting set, we can form four primitive idem- 
potents 

Pi±Pi± = (^)^(1±F,)(1±F,) (200) 

(201) 

We shall denote these by 

Pa^Pi±Pj± (a = 1,2, 3, 4) (202) 

where 

a = 1 — 

a — 2 — > i+.j^ 

a = 3 — ^ 

a = A — > i-J- 
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The triplets of nontrivial commuting F's are (the numbers refer to the index of the 
matrices, see Appendix 0): 

6,9,16 6,4,15 6,5,14 9,2,13 9,3,12 
16,8,11 16,7,10 15,2,11 15,3,7 10,4,12 
2,10,14 11,5,12 14,3,8 13,5,7 8,4,13 

For building P^s we can take any pair from a given triplet, and achieve the same result 
(get the same PaS for the triplet); this is because 

Pa = (^f(i±r,±r,±r,r,) (203) 

and 

r^r, = ±Tk (204) 

Furthermore, each F appears only in three of the above triplets, and we can transform 
from one triplet to another in this "trio" of triplets, by a unitary transformation, thus 
transforming between different sets of Pa's. If we denote the trio of triplets in which a 
F matrix (say, "a") appears as 



b, c 
d, e 
f,9 



then we find that we always have a Pauli type relation of anticommutators among the 
remaining pairs, for example such as: 

{b,d} = {bj} = {dj} = 
{6,e} = {b,g} = {e,g} = 

here we took "b" as the "pivot" (i.e. "b" is used for Tb). This enables us to take T{, 
and transform the second triplet into the third: a,d,e ^ a, f, g (it can be done equally 
well with "c" as a "pivot"), thus getting a transformed set of Pq's. This way we can 
take a triplet in a trio, and use it to transform between the other two triplets. Looking 
at the existing triplets we see that we can transform from any triplet to any other by 
a suitable choice of unitary transformations, i.e., all sets of P^'s are equivalent up to a 
unitary transformation. 



39 



A.2 Nilpotents 

We can produce non-primitive level nilpotents by taking 



n, 



^3 



4 = ±i{r,,r,.} (205) 

and requiring that {Fj, Fj} = 0. 

Now we turn to primitive level nilpotents. In a certain representation we have — 
^kk: {k = 1, 2, 3, 4) where e^fe is the matrix which has a "1" in the A;*'* row and k*^ column, 
and all other places are zero. Since the matrices e^^, (n ^ m) are nilpotent, (we have 
twelve of these), we get for any Dirac matrix A: 

enAekk = ocikCik (206) 

where ttj^ is the element on the i*^ row, k^^ column, in A. Once we have these twelve 
nilpotents, others can be built from them. If we want to represent the nilpotents with 
the help of the F matrices we need to have 

^ikj Pi^kPj 

Nl^ = ArfeP.AFfeP,- = {i^j) (207) 

this is because PjPi = always. We must guard against Ni^j being zero, therefore we 
must not have [Pi,Ffc] = or [Pj,Ffc] = and i ^ j. In this case, for i = j, the 
representation is diagonal. Actually if we take a closer look at Pj, Pj, for example. 

Pi = Pn±Pm± Pj = Pn±PmT (208) 

then if {Tk, F^} = 0, FfcP^± = Pm^^^k, and P^ flips sign. We can summarize the 
nonzero useful nilpotents as follows: 

1. 

{Ffe,F„} = ; [Ffe,F„] = =^ 

Pi^kPj — Pn±Pm±^kPn±Pm^^ — ^kPn±Pm^^ — ^kPj 
Pi'^kPj' — Pn^:Pm±^kPrm^PmT — ^kPn^^Pm^ = ^kPj' (209) 
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2. 



{rfc,r^} = ; {Tk,Tn} = 



pr P p p VP p VP P VP 

thus showing that there are Nikj 7^ 0. The first entry results in four nilpotents which are 
equivalent to the four generated by the second entry up to a sign. Since both cases are 
equivalent we shall concentrate on the second case. We find that any two commuting F's, 
have exactly four mutual anticommuting F's in common. Taking the two commuting 
F's as a pair from a triplet, then for each pair the four anticommuting F's are distinct. 
Therefore we get a "structure" which can be illustrated by a specific example; for the 
triplet P2,-Pii,-Pi5 the relevant structure is : 



Plbi P2 



Tg, Fio, Fi3, Fi4 
Ts, Te; Ty 

Ts, Ts; ri2, F16 



The P's can be thought of as idempotent generators , and the F's as nilpotent generators. 
The structure covers all F's. From the two P's in a row we can build the four P^'s. We 
may take any F from that specific row to produce the four nilpotents (all F's in a row 
produce the same nilpotents). The three rows produce four distinct nilpotents each, 
summing up to the basic twelve required for 4 x 4 matrices. The important thing to 
remember is that by a unitary transformation we can transform from one triplet to 
another, and thus from one structure to another, hence all structures are equivalent up 
to a unitary transformation. A special feature of the above structure {P2, Pu, P15) is 
that it produces the idempotents and nilpotents characterized by the matrices Cjj. 

A.3 Nilpotents for N± 

Since we are looking for which gives us a positive definite probability density, Eq. 
(pTf) , we need to include an "anchor" in the form of 7°, which assures a positive 
part, and up to one more term (which should be smaller or equal to the anchor since 7 
matrices are orthogonal). Now, considering the form of the currents in Eq. (^), 7° = F2 
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is the only anchor possible because of the appearance of 7° in ip, so we have to look at 
the row in the structure where r2 resides and calculate the nilpotents (all other rows 
will not produce nilpotents with this anchor). We take for example the structure: 

PgjPig '■ Tsj ri2, Tis 
Pq, P9 '■ r?, Ts, TiOj Til 

-Pq) P16 '■ Fs, ri4, 

and focus on the first row. Now we generate nilpotents A^j, (i = 1, 2, 3, 4) 



Ni = r2P6--Pl6+ = 


(^)'(r2 + ^r3 


+ «ri2 


-ri3) 


-^2 — ^2P6-Pl6- — 


(^)'(r2 + ir3 


-iri2 


+ ri3) 


N3 = r2P6+-Pl6+ = 


(^f(r2-.r3 




+ ri3) 


A^4 = r2P6+-Pl6- = 


(^m-ir3 




-ri3) 


In terms of the 7 matrices we get 











(^)W-y + 7V + 7V) 






(^r(7°-y-7V-7V) 






(^)'(7° + y + 7V-7V) 




7V4 = 


(^r(7°+y-7V+7V) 


(211) 



We consider these results as given in a specific reference frame. The Lorentz invariant 
form of these nilpotents is: 

N[ = (7-i(^-)) + (7-/(^+))75 
K = (7-^^^-^)-(7-^^^+^)7^ 
= (7-^^'+^) + (7-^^'-^)7' 
K = (7-^^'+^)-(7-/^'-^)7' (212) 

where we used for the four-vectors Z*^^^-' the notation 

= (/o,-/o,0,0) 

= (/o,/o,0,0) (213) 
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and Iq is some number. 

The requirement for a positive definite probability density implies a reduction in half 
of the participating matrices in each nilpotent combination of Eqs. we can do this 

by reverting to the non-primitive level nilpotents: 

. (7-/(^-)) 

^ (7./(i+))75 

(7.t)(l + 75) 

- -(7-s«)(l-7^) 

. (7-/(^+)) 

. (7-/(^-))7' 
. (7.t)(l-75) 

. -(7-3«)(l + 7') 
where the four- vectors s^^\ and t have been used 

t = (to, 0,0,0) 

= (0,si,0,0) (215) 



AT 


1 AT 


1 

2(7 -7 ) 


A T 


A T 


1 / 5 1 1 5\ 

2(7 7 +7 7 ) 


N, 


+ N3 = 


1(7° + 7V) 


Ni 


-iVg = 


|(-7^+7V) 


Ns 


+ A^4 = 


1(7° + 7^) 


N, 


-^"4 = 


|(7°7'-7V) 


N2 


+ = 


1(7° -7°7') 


N2 


-N^ = 


l(-y-7V) 



(214) 



s 



and tQ,Si are numbers. Going to non-primitive level nilpotents, we get twice as many 
nilpotents (eight), which come in pairs differing in sign. Each pair can be a good 
candidate for N_^., N_. From the nilpotents in Eqs. ( |214D , only (7-/'^^''=)) and (7 ■ t)(l±7^) 
are valid for a positive definite probability density (having (7 ■ t) as a Lorentz covariant 
positive definite "anchor"). Furthermore, it seems that (7 ■ l^^^^) and (7-t)(l ±7^) 
should not be treated equivalently, since they are derived from the same row in the 
same structure. 

Carrying out the procedure of finding nilpotents on the twelve available structures, 
(there are fifteen triplets, but three of them include r2 as an idempotent generator 
instead of a nilpotent generator), we find the following nilpotents. We mention only the 
relevant row in the structure containing and giving rise to nilpotents which have a 
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chance of being positive definite. The generating idempotents are in brackets 



P9, {Pq, -Pie) 


: r2 




ri2,ri3 - 


- (7-i(^^)), 


(7 




Pi5, (-P4, -Pe) 


. T2 


Ts, 


r7,rn - 


- (7-/(^^)), 


(7 


t) ± 


Pu, {P5, Pq) 


: r2 


Ts, 




- (7-^(^^)), 


(7 


t) ± C7is(3) 


P9,{Ps,Pl2) 


: r2 


Te, 




(7-i)±(7isW, 


(7 


i) ± i75 


Pii, {Ps, Pw) 


T2 


r5, 


ri2,ri5 - 


- (7-^(=^^)), 


(7 


t)(l±7^) 


Pio, {P7, Pie) 




r4, 


ri2,ri4 - 


- (7-^^^^)), 


(7 


i)(l±7^) 


Pl5,{P3,P7) 


. r2 


r4, 




(7 • t) ± 


(7 


/(2±)) 


PiOj {P4, P12) 






ri4,ri6 - 


(7-t)±(7is(2), 


(7 


t) ± t-f^ 


Pll,{P5,Pl2) 


: r2 






(7-i)±(7is(3), 


(7 


t) ± i75 


Pu,{P3,Ps) 




r5, 


r6,rio 


- (7-/(^^)), 


(7 


t) ± ats(i) 


Pl3,iP5,P7) 


■ r2 


r4, 




- (7-/(^^)), 


(7 


i) ± 




■ i 2 




i 7, i 9 


17 ' J) 


17 












(0,0,±S2,0) 














(0,0,0, ±53) 














(^0, 0, ±lo, 0) 














(/o,0,0,±g 







B Appendix - Solutions of the equations 
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B.l Solutions of the transverse equation 

We consider plane wave solutions that satisfy the transverse Schrodinger-Stueckelberg 
equation as well, of the form 

4>Lrt{p) e-'Ar+v-u^^p^t) (218) 

where u±{p,t) is a four-spinor dependent on the specific momentum and vector t^. We 
observe that Eqs. (|5^,(p5|) commute with the operator 

h{p,t) = j^icrpt)^^ (219) 



This operator was shown by Horwitz and Arshansky |23| to correspond to helicity in the 
frame where = (1, 0, 0, 0), 

j^^iapth' ^ (220) 

so we shall call h{p,t) the extended helicity operator. An interesting feature of the 
extended helicity operator is that in the reference frame where p^ = {po, 0, 0, 0), we get 

jti^^Pth' ^ (221) 

This is because in this frame | p± \= +pov^- The operator h{p,t) can be decomposed 
into two partial operators h{p,t) = h^{t) ■ hp{p,t) which are 

h,(t) = z(7-t)7' (222) 
hpip,t) = ^(apt)(7-t) (223) 

and where they all commute. The operator h^{t) when looked upon in the frame where 
t^j = (1, 0, 0, 0) is seen to be 

hit) Z7V=(^^ o) (224) 

It exchanges the left and right handed parts of the wave function, and gives them a rela- 
tive phase factor. The left part is rotated counterclockwise by 71/2 in the complex plane 
and the right part is rotated clockwise by the same amount. We call this operator the 
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chiral precedence operator because its eigenvalues indicate which chiral part precedes 
the other relative to a counterclockwise rotation in the complex plane. This is an inter- 
esting feature of the solutions of the equations of motion since it expresses an inherent 
broken symmetry between the left and the right parts of the wave function. Moreover, 
it gives physical meaning to the relative phase of the wave function's components, and 
creates a coherence which is preserved through the evolution. This feature of the wave 
function is independent of its momentum. 

We assumed a four-spinor as the basic structure of the solutions, therefore we can 
use two of these operators to characterize completely the four available solutions. The 
chiral precedence operator (for this purpose the product of helicity operator and chiral 
precedence operator has the same effect) also transforms Eq. (^2]), the equation of mo- 
tion, into Eq. (|60|), the equation of motion with interchanged nilpotents, and vice versa, 
therefore relating ip± to ip± 

t{l ■ t)7V± = ¥>± (225) 

This way, the structure of the main field (and auxiliary field) is set to conform to the 
eigenvectors of the chiral precedence projection operator 

P,3± = ^(l±z(7-t)7') (226) 

where Ph^± = Phs±- This projection operator is connected to the projection part of one 
of the nilpotent options discussed in Appendix Eqs. (1^), and is non-equivalent to 
the nilpotents we use, since we cannot transform from one to the other by a unitary 
transformation. 

The solutions are seen to be 

CIW ) (227) 

where the superscript ± denotes the eigenvalue of the eigenvectors when operated upon 
by the chiral precedence operator. When Ph5+ operates on C±(^) it gives 1, and on CZ(^) 
it gives 0. The opposite happens for P/ig-. In the special frame where = (1,0,0,0) 
the solutions are 

^ ( f-' ) '-^-( I- ) ^'''^ 
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and can be taken as two independent vectors, for example 

Since each spinor can be decomposed and characterized by the hehcity projection 
operator 

P,± = ^(l±^(apt)7^) (230) 

where P^^ = Ph±, and by the chiral precedence operator, we use them to characterize 
the four available solutions for a specific momentum. 

In order to find the general form of the solutions, we consider the projection operators 
Ph± and P/i5± which commute. We can form the following non-normalized solutions 

<(p,t) = i(l + e,^(apt)7^)Cl(t) (231) 

where CI(^) ai'e boosted four-spinors composed of stacked two-spinors as in Eq. ( |228| ), 
and is ±1 depending on the eigenvalues of the extended helicity operator. We define 
the normalization constant N{p, t) so that 

iV2(p,t)<^'(p,t)(7 . t)vri{p,t) = 6rr'6ss' (232) 

We turn now to compute the total integrated probability density, which should be 
unity. 

d^X p±rt{x) = / d'^X {(f)±rt{x){^ ■ t)(f)±rt{x)) (233) 



A wave packet is 

2 

(l>±rt{x) = [ d^pN{p,t)J2c{r,s,p\p-t) u''f{p,t)e-'^^+'P^ (234) 

rs 

where c(r, s,p'^,p ■ t) are the weights of the wave packet's components. Performing the 
d'^x integration in Eq. ( P33|) and one d'^p integration, we obtain 

Jd'xp^rtix) = {2Tif j d^p N\p,t) J2 c*{r\s\p\p-t)c{r,s,p\p-t) 

r,r',s,s' 

x<^'(p,t)(7-t)wr(p,t) (235) 

Using Eq. (p32D , we find (one can say that Eq. ( p36| ) results in a constant depending on 
tfj_ rather than 1, and the integration over dH 5{t^ + 1) is normalized to 1) 

j d'xp^rtix) = {271^1 d^'p J2 I c(r,s,/,p-t) |2= 1 (236) 
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We compute the total integrated currents on all space-time and all t^, from Eq. (pB]), 
by using j^^ = -^{(j)±rt{l ■ t) (9^ (j)±rt)- Then we have 

J±f. = I d^x dH 6it^ + 1) J±r^.ix) = I d^x dH 5{e + 1) ^{4>±rt{x){^ ■ 



(2vr)^ / d'p dH 5{t' + 1)^ E I ^^P'^P- \" 

r,s 

(237) 



so the integrated current is 



J±, = (^) (238) 

Now, p^fj^ is space-like, and so is j±^. Since we perform the integration over t^, thus 
having various components of p^, it is possible to acquire a time-like J_l^ current, as in 
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The equations of motion, Eqs. (|64D , (|65D , describe the evolution of the state of a four 
dimensional momentum and currents, existing in the plane orthogonal to the time-like 
vector t^. Considering the frame where = (1,0,0,0), p±fj_ becomes Pi and we are 
dealing with a space-like current. The appearance of as a four- vector in the theory 
is purely for reasons of keeping the invariance of the theory. Another manifestation of 
this aspect is presented in Section VII, where we find the Hermitian term o't'^ f^^, which 
transforms to a*-' in the frame where = (1, 0, 0, 0), and that is why there is no problem 
of Hermiticity in the second order charged equation. The momentum of an event can 
be thought of as being partitioned by into a longitudinal part and a transverse part. 
Each part evolves according to a different evolution equation, and they refer to the same 
vectors p^,tf,. 

B.2 Solutions of the longitudinal equation 

For the longitudinal equation we consider the plane wave solutions of the form 

(PiUp) e-'^^+'P^u\\{p,t) (239) 



Again the equations of motion, Eqs. (|69D,([70D commute with the helicity operator from 
Eq. ( ^19| ), so it can be used to characterize the solutions. Two more operators that 
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commute with the equations of motion are 

htit) = (7-t) (240) 
KdP^t) ^ (7.t)^(apt)7^ (241) 

where /ip^s = ht{p,t) ■ h{p,t), and they all commute. Notice that 

ht{t) = i-f^h{t) 

The characterization is by the projection of helicity Eq. ( pSOj ) and the projection 
operator 

P,,± = ^(l±(7-t)) (242) 

The operator ht{t) in the frame where = (1,0,0,0) is seen to be just 7° which 
is the parity operator; we denote it as the extended parity operator. Therefore the 
characterization of the solutions of the longitudinal equation is done by extended helicity 
and extended parity. The extended parity operator also transforms Eq. (|66D into its 
counterpart with nilpotents exchanged, and vice versa. We obtain the relation 



(7 ■ t)^|| = y^ii (243) 



The solutions are seen to be 



where the superscript ± denotes the eigenvalue of the eigenvectors when operated upon 
by the extended parity operator. In the special frame where = (1, 0, 0, 0) the solutions 
are 

cr - ( n ^" ^ ( -r ) 

and can be taken as two independent vectors. The general form of the non-normalized 
spinor solutions is 

uY{p,t) = 1(1 + e,^(apt)7^)C|y(t) (246) 

where has the same meaning as in the transverse case, and r denotes the eigenvalues 
of extended parity. 
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As for the transverse case we set the normahzation condition on the probabihty 
density 

J d^X p\\rt{,x) = J d^X (0||^t(x)(7 ■ t)(t)\\rt{x)) (247) 

so that 

N\p, t)Mf (J9, t)(7 ■ t)uY{p, t) = Srr'Sss' (248) 

A wave packet is 

0||.i(x) = / d^pN{p,t)Y,d{r,s,p\p-t) (p,t)e-^w-+^P- (249) 

r,s 

and we get after the integration over d'^x and one over d'^p 

jd^xp\U^) = {2t,Y j d'^p N\p,t) d*{r\s\p\p-t)d{r,sy,p-t) 

r,r',s,s' 

xu]^''{p,t){^ ■t)uY{p,t) (250) 

Using Eq. (|2l8D we find 

/ d'xpw^tix) = {2nY f d'p Y: I d{r,s,p',p-t) \'= 1 (251) 

As for the integrated currents using, ji^^ = —21^(011 (7 ' ^) we find 



J||^ = J dS dH 6{t' + 1) J\\r^{x) = J d'x dH 6{t' + 1) _(0||,,(x)(7 ■ t)9||^0,|,i(x)) 

= (27:Y j d'pdH6{e + lf-^ Y \d{r,s,p',p-t)\' 

(252) 

so the integration of the current (over as well) is 

■^11. = (^) (253) 

In a frame where = (1,0,0,0), the zeroth component of the currents is the only 
one which survives due to the vanishing of d±i in this frame, so the currents are time- 
like. In accordance with Lorentz invariance this is true in any reference frame. As in 
the transverse case, J||^ may be time-like, depending on the relative weights of the wave 
packet's components. 
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C Appendix - Discrete symmetries 



When dealing with the discrete symmetries of the theory it is necessary to know how the 
pre-Maxwell field a^t, (a = 0, 1, 2, 3, r) transforms. We require that the free kinetic term 
of the pre-Maxwell field (see |^^), ^Xf^'^fa/s, be invariant. Furthermore, we define some 



new forms of solutions for both the transverse and longitudinal equations, which are just 
the ones defined in Appendix but transformed, and are helpful for a full discussion. 
The basic form of the solutions of the transverse equations are given in Eq. (|227|) ; we 



redefine them by choosing ^ to be an eigenvector of —ia with eigenvalues of ±i 

(254) 

Multiplying each component of C± (^) by —iaH^, and each component of CZ(^) by —ia^t^, 
and exchanging by , we get 

Vtit) ^ ( ) ',!(«) - ( ) (255) 



The basic form of the solutions of the longitudinal equations is given in Eq. ( 244 ), by 
multiplying both components of C|f (^) by (x^'t^, and C|~(^) by — cr'^t^, and exchanging ^"^ 
by we get 

We now proceed to the treatment of parity, charge conjugation, and r reversal. 
C.l Parity 

The transverse charged equations of motion exhibit a symmetry under the inversion of 
parity. The transformation is 

X —X 

P ^ -P 

t ^ -t 

a(x) a(x') = -a(a;) (257) 
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Applying it to the transverse equations of motion we get 

(258) 

Multiplying by 7° from the left we get back the original form of equations 
(am)7Vi,,^,_,(x') = s+M(7-t)70xi.to,-t(^') 

(am)70xU,-t(^') = «-(7-^)(^9. + ea07Vx.*o,-t(^') (259) 
and the wave functions transform (up to a phase) as 

'Plrti.^) = 7V±rto,-t(a;°, -X) 

x!rt(^)=7Vrto,-t(a:°,-x) (260) 

where the superscript P denotes parity. The same is true for the charged longitudinal 
equations of motion, where we find the wave functions transforming as 

(t>\\rt{^) = 7V||rto -t(2;°, -X) 

XH(^)=7°X||rto,-t(^°,-x) (261) 

The solutions of the transverse equations are constructed from a helicity projection 
operator, and a C±(^) part. A similar situation exists for the longitudinal equations. We 
check both components. 

When we check the action of parity on the C± (^) parts of the solutions of the trans- 
verse equation, we find 

PCl(t) = 7°Cl(io,-t) = Clit) 

PCl(t)^7°Cl(to,-t) = a{t) (262) 



and 



P7y+(i) = 7V(^o,-t) = r±{t) 

P77l(i) = 70771(^0, -t) = vUt) (263) 
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On the other hand, for the solutions of the longitudinal equations, we find 

PC|fW^7°C|f(^o,-t) = vtit) 

Pq{t)^^\{to,-t) = v^it) (264) 

and 

pq{t)^^\{to,-t) = C|f(t) 

P^{t) = j\{to,-t) - q{t) (265) 

We can make a combination of states for which this transformation makes them trans- 
from one into the other. This combination is useful also for charge conjugation and r 
reversal. The combination is 

P{q{t)±vtit)) = (CifWi^ifW) (266) 
Doing the same thing for the transverse case one obtaines 

p{Ct{t)±vl{t)) = {Cl{t)±vt{t)) 

P{Cl{t)±vt{t)) = {Ct{t)±Vl{t)) (267) 

On the other hand, when we check the action of normal parity on the helicity projection 
operator part of the solutions, we find 

7°^/.±(Po, -P, to, -t) = PhTiP: t)7° (268) 

therefore the hehcity flips sign. 

Another type of parity is what we call generalized parity, denoted by V, including the 
regular parity and a time inversion. Usually the discrete symmetry under time inversion 
is treated separately and differently from that of parity, as an anti-hermitian operator. 
In our theory time is just another dimension distinguished from the space dimensions 
by the metric tensor, and the evolution is governed by r. Therefore the role of time 
as described by Wigner is transferred to r, the time symmetry becoming much simpler. 
This generalized parity transformation is 
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a^, -a^ (269) 

When we come to define the way the wave functions transform we have a few options. 
The charged form of Eqs. (|6^ ) , (|65|) 



(7-t)(am)<^x = s+Mx± 

{-f ■ t){am)x± = 2s.{idr + ear)^± (270) 

permits two compensating operators for the transformation of Eqs. (|269| ), 7^, and (7 ■ t). 
We prefer to use 7^ as a generahzed parity operator because it is not coupled to any 
four-vector. Applying the transformation we find 

-(7-t)(ant)0l,,_,(x') = s+Mxl.,_i(x') 

-{j-t){am)x'^^_,{x') = 2s_{zdr + ear)<P'^^^^t{x') (271) 
and after multiplying by 7^ from the left we obtain the usual form 
(7-t)(am)7^0'^,^„,(x') = s+M^'x'±r,-ti^') 

{^■t){amh^x'±r,-ti^') = 2s^{tdr + ea,h^<p'^,^^,{x') (272) 
In this case, the wave functions transform (up to a phase) as 

(Plrtix) = 7V±r,-t(-a;) 

Xlrti^) = l'X±r,-t{-x) (273) 

A similar transformation holds for the longitudinal equations of motion. As can be seen 
by operating on the solutions of transverse and longitudinal equations, we get 

m,i|W = 7^Cl||(-t) = +Cl||(t) 

^C7||(t)=7'Cl,||(-t) = -C7|,(t) (274) 



and 



Pr/+||(t)=7V,i|(-t) = -<||(t) 

T^Vliiit) ^ l'vl,iii-t) = +Vui{t) (275) 
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Therefore, the generahzed parity operation on the transverse equations, brings the set 
of solutions into themselves. 

Here too, we can use combinations which behave under generalized parity the same 
as shown for the regular parity, but now the behavior of the transverse and longitudinal 
solutions is the same. For the transverse and longitudinal cases we have 

nCl||(t)±r/I,|,(t)) = +(Cl||W±r/I,||(t)) 

nCl,fi)±vUt)) = -(Cl,||W±<||(t)) (276) 
Note that the combinations 

C(t)-vHt) (277) 
are pure right handed, and the combinations 

C-(t)+r?+(t) (278) 

are pure left handed, for both ± and ||. The extended helicity operator remains un- 
changed under V, so there is no helicity flip in this case. 



C.2 Charge conjugation 

The theory is symmetric under charge conjugation. We define the charge conjugation 
operation, denoted by C, as the transformation necessary to bring the charged equations 
of motion to the same form only with the sign of e reversed. First conjugate the equations 

{l-tna^''*{-P,-eaM<f>'Ux') = s+MxlM) 

(ry.t)*(a'^'^*(-P^-ea^)Qx'lM) = 2s_(-i9, + ea,)0'_:,,(x') (279) 
then make the substitution 



—ttj 

—T 



(280) 
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which is required to bring the equations to the correct form. This is consistent with the 
invariance of the free kinetic term of the pre-Maxwell field IA/^'^/q/?. The next step is 
to multiply them by («7^7^) = (7^(77°) where C = ^7^7°, because 

iJ^'J^{a^'^*){i-f^J^y^ = -a^"" (281) 

One then obtains 

(7-t)(am)z7S'0'j:-.t(a:') = s+Mt^'^\'l_^,{x') 

(7-t)(am)27W:-.t(a;') = 2s4tdr-ea^)t^'^^'*_^,{x') (282) 

where 11^ = + ea^. Taking —a^ brings us back to the original form. Therefore, 

the wave functions transform (up to a phase) as 

<Pirti^) = «7'7V*±,-r,t(x) 

Xirti^) = ^l'l'xl,-r,t{^) (283) 

The same persists for the longitudinal equations. This transformation includes a r 
reversal for its consistency. When the pre-Maxwell field radiation equations [^| are 
taken into consideration, 

d,r^ = ef = ep (284) 

a symmetry breakdown is observed in the second equation due the need to change the 
sign of the field. This is a consequence of the Schrodinger type equation used. Here 
p is a scalar. One could think of a Klein-Gordon generalization for which jV oc ilj*drip 
for which this problem would not occur (note that the same question arises in the non- 
relativistic Schrodinger quantum theory in interaction with the standard Maxwell field). 

To see the consequences of performing charge conjugation to the solutions of the 
transverse equation, we define solutions as in Eq. ( p31| ), using rj\_{t), 

w^lip^t) = ^(1 + e,j^^{apth')vl{t) (285) 
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The subscript s in Eq. ( p85| ) denotes the eigenvalue of the extended hehcity operator. 
Considering a free event, we have for the hehcity projection operator 

ZT^'^/ti = PhT^l'l^ (286) 

and since 

- za'e- = -r 

-ta'r* = +t (287) 
for the Cl(^) P^-rt of the solutions, we have 

cent) = -viit) 

cent) = +Vlit) (288) 



and 



Cvtit) = +C^(^) 

Cvlit) = -CZ{t) (289) 



so we can use the combination 



CiClit)±v-At)) = -{Cl{t)±vt{t)) 

C(Cl(t)±r/l(t)) = HClit)±Vl{t)) (290) 

These combinations have mixed eigenvalues of the chiral precedence operator, positive 
for C^{t), f]~l{t), and negative for C]^{t), rj^^it). Therefore for the combination we obtain 

2 2 

Ce-'^^+'P%uf{p, t) ± wn\p, t)) = -ere-'^^-'P%wn\p, t) ± ^^'"^(p, t)) (291) 

where is ±1 depending on value of r, either + or — . The state transforms into 
the other helicity state. Since no space inversion was performed, we get a helicity flip. 
Therefore, the charge conjugated wave function describes an event with opposite charge 
and spin, moving in the opposite direction. As mentioned in Appendix |U.1| dealing with 



parity, these combinations are pure left or right handed spinors. For example, the charge 
conjugated left handed spinor is a right handed spinor with opposite charge moving in 
the opposite direction. 

Repeating these operations for the solutions of the longitudinal equations, we find 
similar results, we have only to change the subscript ± to ||. 
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C.3 r reversal 



Here we discuss r reversal in the sense introduced by Wigner; reversing the evolution 
parameter we reverse the motion of the event, thus creating a current in the opposite 
direction. The operations needed for the r reversal symmetry, denoted as T, result in 
the same operations as those for the charge conjugation as can be seen by performing 
the operation 

r ^ -r (292) 
on the equations of motion. To compensate we have to take 

a^, -a^ (293) 

and applying this to Eqs. ( p7CI| ) while taking the conjugate of the equations, we get the 
transformation properties of the wave functions 

<Plrt{x) = 7'7V*±,-r,t(x) 

Xlrtix) = I'l'xl-rA^) (294) 

As usual, the longitudinal equations have the same transformation properties. We can 
repeat the same arguments given for the charge conjugation. Although the transforma- 
tion looks the same as the charge conjugation one, the reason for doing it stems from 
Wigner's idea of time reversal. However, we obtain a charge conjugated event. 

C.4 VCT invar iance 

The most striking thing about these discrete transformations is the identity of the charge 
conjugation transformation, and r reversal transformation in the sense of Wigner. Con- 
ceptually this is understood using Feynman diagrams. We know from Dirac's theory 
that a charge conjugated electron, i.e. a positron, moves backward in space-time with 
opposite charge. This is exactly what a r reversal means. Seen in the four dimensional 
world (a projection of the trajectory on space-time at a specific r), a r reversal changes 
the charge, direction, and helicity. Concluding that charge conjugation and r rever- 
sal have the same consequences, a full VCT transformation is just a generalized parity 
transformation. We get an event with opposite motion in space-time relative to the state 
before conjugation, with the same helicity and charge. The symmetry under generalized 
parity is a direct consequence of the manifest covariance built into the theory. 
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D Appendix - Dirac matrices 



Our orthogonal matrix basis which satisfies the Dirac Chfford algebra is |^ 



ri = i 12 = 7° r3 = ^7' r4 = z7' 
r, = t^^ r6 = 7V r7 = 7V r8 = 7V 

Tg = 27 V Tio = i7^7^ Til = 27^72 = 7I7273 

ri3 = ^7V7' ri4 = ^7°7V ri5 = 27°7V ri6 = z7V7V 

We use the following chiral representation of the Dirac matrices 



7-1 ° M y= ° 

' I 1 i I -a* 



and 



7^ 



-1 
1 
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